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Abstract 

We consider limits of equilibrium distributions as temperature ap- 
proaches zero, for systems of infinitely many particles, and character- 
ize the support of the limiting distributions. Such results are known 
for particles with positions on a fixed lattice; we extend these results 
to systems of particles on M. n , with restrictions on the interaction. 



1 Introduction 

Consider a physical system consisting of a large number of interacting molecules 
in thermal equilibrium. Equilibrium statistical mechanics accurately models 
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such systems. The fundamental qualitative feature of having fluid and solid 
phases, the latter appearing at low temperature and/or high pressure, can 
even be usefully modelled with classical statistical mechanics. Yet although 
this phase structure has been amply supported by computer simulations [FS] 
there is as yet not a single model, of particles moving in space and interacting 
through reasonable short range forces, in which such fundamental features 
can be proven [RalJ. It is not difficult to model a solid if one uses less physical 
interactions, as in the Einstein or Debye models, but this has not yet been 
achieved with a satisfactory short range model, which could also describe a 
fluid phase. (See however |Rult ILMt IRa2"] .) This is one of the main unsolved 
problems in condensed matter physics \Br\ lUhj . 

One of the difficulties is that to unambiguously characterize or distinguish 
a phase one must use the thermodynamic limit, or, equivalently, uniformly 
control the behavior of the system as the system size grows indefinitely. 
Now consider the grand canonical ensemble for a system of finite size, for 
which the unnormalized probability density of particle configuration u is 
exp —(3[E(u) — \\u\], where E(oS) is the energy of u, \u\ is its particle number, 
A is the chemical potential and is the inverse temperature. If one takes 
the limit — > oo, for fixed A, one easily sees that the probability tends 
to concentrate on configurations u which minimize E(u) — X\u\ and which 
are called ground state configurations. This can be useful, as one can then 
understand the state at large (3 as a perturbation of the energy ground state, 
and try to understand the solid phase from this. But, as we just noted, one 
must be careful with the order of limits; one needs to see the approximation 
of positive temperature states by zero temperature states uniformly in the 
size of the system, and not only is this not obvious, it can actually fail, as 
we show below. 

In other words, a major difficulty in solving this old problem, of success- 
fully modelling the origin of the solid state in terms of short range forces, 
is to control the approximation of low temperature states by those at zero 
temperature (energy ground states) as the system size grows. One way to sys- 
tematize this is to make sense of infinite volume limit "Gibbs states" fip t \ for 
finite (3, introduce the "ground states" //oo,a as the limit points lim^oo fJ,p t \, 
and characterize the "ground state configurations" uj in the support of /^oo,a, 
the smallest closed set of configurations of probability 1. This is the path we 
will take as it makes the control of the limit as f3 — > oo a bit easier, having 
already taken the limit in the size of the system. So /i^a and /ioo,A will be 
probability distributions on configurations u of particles in unbounded phys- 
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ical space. From the above analysis in finite volume we expect /ioo,A to be 
supported by configurations u which in some sense minimize E(u) — X\u\. 
Of course for an infinite system E(u) — X\u\ is typically going to be infi- 
nite, so one must adjust appropriately both the defining characteristic of the 
equilibrium distribution fjLp t \ and the optimization approached as — > oo. 
For np t \ this was solved rather generally many years ago, and we use the 
result below. But for /x^a this was only solved in the simpler situation of 
models with particles living in a discrete space, typically a lattice |Ru2| ISc]. 
Restricting oneself to lattice models is a weakness however if one is trying to 
model the fluid/solid transition, or more specifically to model a solid. Re- 
moving this obstacle is the main motivation for this work. We show that 
under certain reasonable conditions the optimization characterization of the 
support of /ioo, a is the same for models of particles in space as for models of 
particles on a lattice. The proof is somewhat harder, and for good reason: 
it is known that without extra assumptions on the interaction the particle 
density can be unbounded at any positive temperature, a phenomenon not 
possible in typical lattice models. 

Our arguments are necessarily technical since we are forced to deal care- 
fully with limits, but this is justified by the direct relevance of our results 
to matters of importance to physical theory. We note that ground states 
are also used significantly in optimization schemes outside physics; see for 
instance [HRJ. 

2 Convergence to Ground State Configura- 
tions 

First some notation and assumptions. We assume a two-body interaction 
potential U (s, t) dependent only on the separation of the point particles at 
positions s, t in R n , including a hard core at separation 1 and diverging as 
the separation decreases to 1: 



We assume U has finite range R > 1, and that U > — m, for some m > 0. 
Denote the chemical potential by A. 

We denote by Q the set of all finite or countably infinite configurations 
w C R" of particles which are separated by a distance at least 1. By Uj we 




for \s — t\ < 1, 
for \s — t\ \ 1. 
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denote the positions of the particles in u, and by bi(u) the set of balls bi(u}j) 
of diameter 1 centered at positions uuj. For A C W l we denote by Q A the set 
of configurations u = ooa = oo n A, which have all their particles in A. The 
number of particles in uia will be denoted by \u>a\ ■ With the usual topology 
Q is compact. 

Now we will introduce the notion of the Gibbs state. We first introduce 
for every bounded A and u> G Qa the energy: 

H(u) = s ^U(u i ,u j ) + \\u\. (1) 

For two collections of particles, u/ G Qa, w" G f2, we define the interaction 
between them as 

H(u J ',u J ") = J2m,^)- (2) 

and the sum 

H(u/\u") = H{u') + H(u', u"). (3) 

We say that the probability measure /z^a on Q is a Gibbs measure, corre- 
sponding to the interaction U, inverse temperature (5 and chemical potential 
A, if for any finite A and any function / on Qa we have 

/N d/j,^ x (uj) (4) 



f (u A )exp{-[3H(uj A \ujAc)} d-K A (u A ) 



ZA (A,u A c)dnp, x (w). 



Here 

• A c = R n \ A, 

• 7i a is the Poisson measure on Qa, which on the /c-particle subset of Qa 
is just the Lebesgue measure, normalized by the factor e - '^'^-, 

• the partition function 

Z(3 jX (A,ujac) = J exp{-pH(u A \u}Ao)} dir A (uJA) ■ 
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It is easy to see that any such measure gives probability 1 to the set of 
configurations in which no two particles are at distance 1. 

The probability distribution qA,u A c = Qp,\,A,u A c on ^a, given by the den- 
sity Zp\ (A, loa*) exp {— /3H(uja\wa c )} with respect to the measure ha, is 
called the conditional Gibbs distribution, corresponding to the boundary con- 
dition uja c - The equation (TJ]) is called the Dobrushin-Lanford-Ruelle (DLR) 
conditions; see |D1[ ID2[ ID3[ ILRj . (The conditions are a way to replace, for 
systems of infinite size, the usual formula which one uses for finite systems.) 
Any measure obtainable from such a Gibbs state by the limit /3 — > oo will be 
called a ground state. 

Let the set G of ground state configurations be defined as: 

G = {uj G ft : for every bounded Ad" and every uj = (cu' A , uj\c) , 
H(J k \ujkc) - H (w a |waO > 0}. 

This set is nonempty, see |Ra2] . 
Our main result is the following 

Theorem 1 Let fi^ be any limit point of the family of Gibbs states fip as 
— > oo, i.e. a ground state. Then fi^ (G) = 1. 

Theorem [1] holds in more a general situation, when the interaction has 
no hard core, but possess instead the superstability property. (We discuss 
superstability in the last section.) The proof is more complicated, and we 
will not present it here. 

Theorem [TJ is equivalent to 

Theorem 2 Assume that u G G c . Then there exists an open neighborhood 
W of u such that J w dfi^ (a) — > as (3 oo. 

Proof of Theorem [2], Before giving the formal proof we present its 
plan. If uj G G c then the following holds: there exists a finite volume B, inside 
which the configuration uj = (ub,ujb^) can be modified into Q = (u! B ,ub^) 
in such a way that 

A (co) = H (ub\ub c ) — H (ujb\ojb c ) > 0. (5) 
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We will be done if we can find open neighborhoods W, W of the configurations 
co and Co, such that 

lie (W) 

as (3 — > oo. So we need to find an upper bound for pp (W) and a lower bound 
for fjLp (W) . To do this we will use the following simple 

Lemma 3 For every value of the chemical potential A there exists a distance 
p (A) > 1 such that the following holds for all p in the interval (l,p(A)): 

Let M C R™ be any bounded volume and £ G VLm c - any "boundary 
condition". Denote by Qm, p (0 C fin the subset 

a G Qm '■ t wo particles of a are separated by < p, or a particle of a 
is at distance < p from a particle of £ 

(7) 

Then the conditional Gibbs probability qp,M£ (&M,p (0) 9 oes to as (3 oo. 
This convergence, of course, is not uniform in M, but for every M it is 
uniform in £. In other words, 

qp,M£ (fiif \ ^m, p (£)) = 1 - 7 (A M, £, p) , (8) 
where for every M, p the function 7 M, £, p) — >■ as /3 — > 00, uniformly in 

c 

T/ie same statement holds for the subset 

&M,p = {o~ G f^M : ^0 particles of a are at distance < p] , (9) 
since /or ewer?/ £ we /lave fiM,p C ^M,p (0 ■ 

Without the hard core condition Lemma [3] does not hold, and has to be 
replaced by a weaker statement. Our proof of Lemma [3] uses the divergence 
of the repulsion near the hard core. 

The proof of Theorem [2] proceeds now as follows. Let B be the open 
R- neighborhood of B in M n . Due to the Lemma it is enough to consider the 
case of 00 such that u§ ^ Q§, P (\)- 

By an r-perturbation of a finite configuration w G Q we will mean any 
finite configuration x with the same number of particles, such that for every 
particle Wj G w the intersection >cf\b\{wf) consists of precisely one particle 

G x, and dist(zUj, >cf) < r. 
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Now we define the open neighborhood W of u by putting 



(10) 



where Q r (u, is the set of all those r-perturbations x of cu B which also 
belong to fl B . It is immediate to see that if r < p(A)/2 then for every (x, £) e 

|#(x B ) -H(u B )\ <Cr, 
\H (x b , xg\ B ) -H(u b , vb\b) I < C r > 



for some C = C (B) . Let r be so small that Cr < 



AH 
10 ' 



Then, by DLR, 



dfx (x, f ) 



L r (c^ xs ) exp (x B |x As ) } d7r s (x s ) 



f7 r cj,B 



< exp <J -f3 
where 



H ' (u B \u) B c) - 



Z B 
A( W ) 



■d/x (x,0 
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n r (w,B) Z B (x b \b) 



d/i (x,0 , 



fi r (w, xg\ B ) = {xGO r (w,B): x BVB = x BVB } , 

and Z B (xm B )'s are the partition functions. 

In the same way, and recalling the meaning of Co, we put 



(12) 



Without loss of generality we can assume that for the same C and every 

(x,0eiy 

\H (x B )-H(u B )\ <Cr, 
\H (x b , Xb\ B ) ~ H (uj b , ub\b) \ < Cr - 
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Then 



d/i 



r\W,X B \B 



exp {-f3H (xb\xb\b) } dn B (x B ) 



a r [ui,B] 



z v \Xb\b) 



dn (x, 



> exp < —ft 



H(u b \oj B c) + 



AM 



■In 



d7i B (x B ) 
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But the integral f Q \ dn B (x B ) is just the Poisson measure of the 

set Q r (u) , so it is a positive number (not depending on (3). The comparison 
of the last two estimates proves our theorem. □ 

Proof of Lemma [31 Using the above ideas it is straightforward. Let 
% (n) be the maximal number of particles with which any given particle can 
interact. Suppose a particle w\ is p-close to zu 2 . Due to the divergence of the 
repulsion near the hard core, the energy of the interaction of particles w\ and 
ZU2 diverges as p \ 1, so we can assume that U (zui, w<i) > \+i {n) m+1 once 
p — 1 is small enough. But then if we erase the particle zn\, we gain at least 
one unit of energy. The rest of the argument follows the same line as above. 
□ 



3 Counterexamples 

In this section we will explain that some results which one might expect to 
obtain in this area in fact do not hold. 

Let U be some pair potential, which is translation and rotation invariant, 
i.e. U (s,t) = U (\s — t\) . We suppose U to be superstable. Superstability is 
a property which means that the repulsion part of the interaction dominates 
the attraction part, see |Ru3] for more details. The Lennard- Jones potential 
is an example of such an interaction. Let A be some chemical potential. 

Our initial modest goal was to prove that for all reasonable interactions 
U the following holds: 

Statement. Let (3 n — > oo be a sequence of inverse temperatures, going to 
infinity, and let p n be a weakly converging sequence of Gibbs states, corre- 
sponding to the interaction U, chemical potential A and inverse temperatures 
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P n , i.e. p n G Q (U, A, (3 n ) . Then the limiting state floo is supported by the set 
G of ground state configurations. 

For every U, A there exist a pair of constants R<R, such that for every 
uj = {uj k } e G 

inf \uji — ujA > R, (13) 

i<j 

and every ball B r C R n of radius r > R contains at least one particle toi G u. 

If true, these properties would be a zero-level approximation to the or- 
dered structure which is expected (in some sense) to be formed by solids. 
(For the use of these Delone properties in modelling the ground states of 
quasicrystals see [BHZl ED Hi 03 |Ri2l [RWl NGW] .) 

We have to say that we do not expect the above picture to hold without 
extra assumptions, though these assumptions are expected to be mild and 
physically natural. The following results point to difficulties that must be 
overcome. 

Let /i G Q (U, A, (3) be a Gibbs field with inverse temperature (3, and the 
superstable interaction has attractive part. Denote by p M (x) the expected 
number of particles of the field /i in the unit ball centered at the point x G M n . 

Proposition 4 For every A, and U without hard core there exists a state 
G Q (U, A, 0) , such that the function (•) is unbounded on W 1 . 

That statement means that the relation (TlBl cannot hold in general. 

Proposition 5 Suppose that the state jl G Q (U, A, 0) has the density func- 
tion Pp,{-), which is polynomially bounded, i.e. there exists a polynomial 
P (•) , such that pp. (x) < P (x) , x G M n . Then there exists a constant 
C = C(U,X), such that PfL (x) < C. 

The proof of Proposition [5] can be obtained by the application of the 
technique of compact functions, developed by R. Dobrushin in |D4j . see also 
[DPJ. Being proven, Proposition \5\ can be used to deduce the existence of 
the constants R, R above, under condition that the random fields we are 
dealing with have their density functions polynomially bounded (and hence, 
uniformly bounded). 

Proposition H] can be derived from Proposition [5] and the following con- 
struction. We will consider the ID case; the generalization to higher dimen- 
sions is obvious. Let us suppose that U (r) > for r < rj, U (r) < for 
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ri < r < r 2 , and r*i < 1/3, r 2 > 2. Let /„ be the unit segment centered 
at the integer point n G R 1 . Let w^ 1 , to 1 be two configurations in the seg- 
ments 7_i, 7i, and consider the conditional Gibbs distribution q (u; |a7 _1 , zu l ) 
in I , given configuration tu^U vj 1 outside. Let K > be fixed. Clearly, 
there exists a number N (1) , such that if | -ccj 1 1 > N (1) , 1 > iV (1) , then 
27T > E (|u; |) > K. Indeed, there will be a part of the segment Iq, where the 
potential defined by the particles iu _1 U w 1 will be attractive (=negative), 
and the more particles we will have in a7 _1 U w , the deeper this well will be. 
Let now w~ 2 , w 2 be two configurations in the segments i_2, h, and consider 
the conditional Gibbs distribution q (u\m~ 2 , w 2 ) in J_i U Io U I\, given con- 
figuration m~ 2 U w 2 outside. Clearly, there exists a number N (2) , such that 
if \w- 2 \ > N (2) , \w 2 \ > N (2) , then E ((w 1 !) > JV(1) , E ( 1 1 ) > N (1) , 
and so again E (|o; |) > K. Here we denote by u k the restriction of u; on the 
segment Ik- If the number N (2) is not too big, then we have in addition that 
E (\oj°\) < 2K. We can repeat this construction inductively in n. As a result, 
by taking a limit point we get an infinite volume Gibbs state on R 1 , such 
that E (|o; |) > K. If K is chosen large enough: K > C (U, A) , then the so 
constructed state has the function p(-) unbounded, due to Proposition 

In the previous section we presented a proof of our Statement, restricted 
to the case of interactions with hard core. This assumption plays a technical 
role, and with an extra effort it can be removed. 
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